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The phenomena customly described with the standard ACDM model are broadly re- 
produced by an extremely simple model in TeVeS, Bekenstein's (2004) modification of 
General Relativity motivated by galaxy phenomenology. Our model can account for the 
acceleration of the universe seen in SNela distances without a cosmological constant, and 
the accelerations seen rotation curves of nearby spiral galaxies and gravitational lensing 
of high-redshift elliptical galaxies without cold dark matter. The model is consistent with 
BBN and the neutrino mass between 0.05eV to 2eV. The TeVeS scalar field is shown to 
play the effective dural roles of Dark Matter and Dark Energy with the amplitudes of the 
effects controled by a /i-function of the scalar field, called the fi— essence here. We also 
discuss outliers to the theory's predictions on multi-imaged galaxy lenses and outliers on 
sub-galaxy scale. 
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1. Introduction 

As with the start of the last century, a rethinking of fundamental physics is forced 
upon us by a set of experimental surprises, only difference with this time is that 
the whole universe is the laboratory. Einstein's General Relativity together with the 
ordinary matter described by the standard model of particle physics is well-tested 
in the solar system, but fails miserably in accounting for astronomical observations 
from just the edge of solar system to cosmological distance, e.g., fast-rotating galax- 
ies like ours would have been escaped the shallow gravitational potentials of their 
luminous constituents (stars and gas). Standard physics also cannot fully explain 
the cosmological observations of the cosmic acceleration seen in Supernovae type 
la data and the angular scales seen in the anisotropy spectrum of the Cosmic Mi- 
crowave Background Radiation(CMBR) (Spergel et al. 2006). The remedy is usually 
introducing two exotic components to dominate the matter-energy budget of the 
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Universe with a split of about 25% : 74% to the universe energy budget: Dark Mat- 
ter (DM) as a colissionless and pressureless fluid described by perhaps the SUSY 
physics, and Dark Energy (DE) as a negative pressure and nearly homogeneous field 
described by unknown physics. 



1.1. Challenges to dark matter and dark energy 

In spite of the success of this concordance model the nature of dark matter and 
dark energy is one of the greatest mysteries of modern cosmology. For example, it 
has long been noted that on galaxy scales dark matter and baryonic matter (stars 
plus gas) have a remarkable correlation, and respects a mysterious acceleration scale 
ao (Milgrom 1983, McGaugh 2005). The Newtonian gravity of the baryons g fc and 
the dark matter gravity gDM are correlated through an empirical relation (Zhao & 
Famaey 2006, Angus et al. 2006, Famaey et al. 2006) such that the light-to-dark 
ratio 

_9±_ = 9DM + ag^ ^ = (J) 
9dm a 

where ao is a dividing gravity scale, and < a < 1 is a parameter, experimentally 
determined to fit rotation curves. a Such a tight correlation is difficult to understand 
in a galaxy formation theory where dark matter and baryons interactions enjoy a 
huge degrees of freedom. Equally peculiar is the amplitude of dark energy density A, 
which is of order 10 120 times smaller than its natural scale. It is hard to explain from 
fundamental physics why DE starts to dominate the Universe density only at the 
present epoch, hence marking the present as the turning point for the universe from 
de-acceleration to acceleration. This is related to the fact that A ~ a 2 ,, where ao is 
a characteristic scale of DM. Somehow DE and DM are tuned to shift dominance 

2 

when the DM energy density falls below -^q- These empirical facts should not be 
completely treated as random coincidences of the fundamental parameters of the 
universe. 

Given that the dark sector and its properties are only inferred indirectly from 
the gravitational acceleration of ordinary matter, one wonders if the dark sector 
are not just a sign of our lack of understanding of gravitational physics. Here we 
propose to investigate whether the role of DM and DE could be replaced by the 
scalar field in a metric theory called TeVeS. 

2. TeVeS Framework 

TeVeS is a co- variant theory proposed by Bekenstein which in the weak field limit 
reduces to the phenomcnogically successful but non-covariant MOND theory of 
Bekenstein & Milgrom (1983). The co- variant nature of TeVeS makes it ready to be 
analysized in a general setting. 



''Note a mimics the role of the mass-to-light, hence inherits some of its uncertainty. 
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Just like Einstein's theory, Bekenstein's theory is a metric theory, in fact, it 
has two metrics. The first metric g^ u is minimally coupled to all the matter fields 
in the Universe. We shall call the frame of this metric the "Matter Frame" (MF). 
All geodesies are calculated in terms of this MF metric. For example, in a quasi- 
static system like a galaxy with a weak gravitational field, we can define a physical 
coordinate system (t,x,y,z) such that 

dr 2 = g^dx^dx" = (1 + 2<$>)dt 2 - (1 - 2$) (dx 2 + dy 2 + dz 2 ) . (2) 

Here the potential $ = + where the field <p replaces the usual role of the 
potential of the Dark Matter. 

Another metrics of TeVeS g^ u has its dynamics governed by the Einstein-Hilbcrt 
action S g = 16 * G / d 4 xy^gR, where R is the scalar curvature oig^. We shall call 
the frame of this metric the "Einstein Frame" (EF). It is related to the MF metric 
through g^ = e~ 2( t '{g t _ LU + A^A V ) — e 2 ^ ' A^A V (notations of tildes here are opposite of 
Bekenstein), which involves the unit timelike vector field A^ (can often be expressed 
as (v/— 500) 0: 0, 0) f° r galaxies or FRW cosmology), and a scalar field <j>, which is 
governed by the action S = J d A x^J~^gL, where according to Bekenstein (2004) and 
Skordis ct al. (2006), the Lagrangian density 



C = -A+ ' 



16ttG 



dusk 



(3) 



{p v - A^Ae v ) 4>^, v = ¥ = (4) 



where A is a constant of integration, equivalent to cosmological cosntant, and V is 
a free function of fist, which is an implicit function of the scalar field <p through 

dV 
d^isk 

By picking an expression and parameters for the scalar field Lagrangian density 
C or potential V, one picks out a given TeVeS theory. b Note in all above G = 
(1 — K/2)G® is a to-be-determined bare gravitational constant related to the usual 
experimentally determined value G© w 6.67 x 10~ n through the coupling constant 
K of TeVeS (Skordis, private communication). 



3. Connecting galaxies with cosmology 

Bekenstein's original proposal was to construct the Lagrangian density with C as 
one-to-one function of the nsk (cf. Skordis et al 2006). Such one-to-one construction 
has the drawback that the Lagrangian necessarily has unphysical "gaps" such that a 
sector is reserved for space- like systems (e.g. from dwarf galaxies to the solar system 
in < [isk < Mo) and a disconnected sector is reserved for time- like systems (e.g., 
expanding universe in /isk > 2/Uo)- While viable mathematically, such disconnected 
universe would not permit galaxies to collapse out of the Hubble expansion. The 
particular function that Bekenstein used also result in an interpolation function, to 



b These notations are directly related to Bekenstein's notations by ^i$k = 8717//A;, ¥ = -r|y. 
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be computed from a non-trivial implicit function of the scalar field strength, which 
is found to overpredict observed rotation curve amplitudes when the gravity is of 
order an (Famaey & Binney 2005). 

In an effort to re-connect galaxies with the expanding universe Zhao & Famaey 
(2006) proposed to construct the Lagrangian as a one-to-one function of the scalar 
field (j) through ¥, where ¥ > in galaxies, and ¥ < for cosmic expansion. This 
way allows for a smooth transition from the edge of galaxies where ¥ ~ to the 
Hubble expansion. Zhao & Famaey also suggested to extrapolate the Lagrangian 
for galaxies to predict cosmologies to minimize any fine-tuning in TeVeS. The coun- 
terpart of the ZF model in dark matter language would be Eq. 1 , which they used 
to fit to rotation curves, and found that both a = 1 and a — give reasonable fits, 
with some preference on the former. 

Our aim here is to check whether the suggestions of Zhao & Famaey (2006) lead 
to reasonable galaxy rotation curves and cosmologies. To minimize fine-tuning, we 
consider an extremely simple Lagrangian density governing the scalar field 



J 87rG© a o exp(-0 o ) 

In quasi-static systems \/|¥| = |V</>| exp(— 0), where the constant 4>o is the present 
day cosmological value of the scalar field <j>. With this, the Poisson's equation reduces 
to —V ^ l Vex P(^o~^)l — A-kG^p = — V • gb- So in spherical approximations we 
have 

a{4>) = a exp(<j> - <j> ). (6) 



|V0| a(4>) ' 

Clearly the above Lagrangian or TeVeS Poisson's equation for the scalar field essen- 
tially recovers eq. (1) in the a = case in the Dark Matter language if we identify 
that the scalar field V</> — > gDM, hence playing the role of gravity of the dark matter 
9dm at present day when <\> — fo. Interestingly the characteristic acceleration scale 
a{4>) varies with redshift together with the scalar field <fi(t). 

Tabic 1. Comparison of the pros (+) and cons (-) of LCDM vs. TeVeS in various scales 
Data references 

Pioneer Anamoly - + Sanders (2006) 

Rotation Curves HSB/LSB ++ Zhao & Famaey (2006) 

Lensing by Ellipticals ++ +/- Zhao, Bacon, Taylor, Home (2006) 

Dynamics of X-ray Clusters ++ +/- Sanders (2003), Pointccoutcau & Silk (2006) 

Hubble expansion and CMB ++ +/- Zhao et al. (2006), Skordis ct al. (2005) 



A summary of how well TeVeS /MOND or CDM fits data on all scales is given 
in Table 1. To illustrate, two sample fits to rotation curves of a dwarf galaxy and a 
high surface brightness spirial galaxy are shown in Fig. la, b, including the possible 
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effects of imbedding the galaxies in a large neutrino core. We also repeat the ex- 
cercise of Zhao, Bacon, Taylor & Home (2006), and fitting the lens Einstein radii 
with Hcrnquist models in a a = modified gravity. We show in Figlc,d that the 
CASTLES gravitational lenses (mostly high rcdshift ellipticals) are mostly consis- 
tent with TeVeS predicted Einstein ring size within plausible uncertainties of the 
mass-to-light ratios. Note that the critical gravity (c 2 /Di)(D s /Di s ) is always much 
stronger than 10 _10 m/s 2 at Einstein radii of elliptical galaxy lenses, so the Einstein 
rings are insensitive to MONDian effects, hence insensitive to ao. Some of the out- 
liers are known to in galaxy clusters (RXJ0921 and SDS1004), where a neutrino 
density core of a few times 1O _6 M kpc -3 might help to reduce the discrepency. 
Given that the a = model is reasonable consistent with spiral galaxy rotation 
curve data and Einstein rings of high rcdshift ellipticals, next we wish to study 
cosmology in this TeVeS model. The important thing to note here is that the cos- 
mological constant A is set to zero, so the zero point of the Lagarangian coincide 
with where the scalar field is zero. 



4. Hubble expansion and late time acceleration 

TeVeS is a metric theory, the uniform expanding background can be described by 
the FRW metric. Assume a flat cosmology with a physical time t and scale factor 
a(t), we have 

ds 2 = -dt 2 + a 2 (t) [d X 2 + X 2 (d0 2 + sin 2 Odcf)} . (7) 
The Hubble expansion can be modelled with 

P4> + Pb + Pr = &7rG<sT (8) 

where the first term is the scalar field effective energy density = ¥^ — C = 

^pexp(50) [jft^j (8-7rG©a exp(— O )) 1 in the matter frame. The correction fac- 
tor 

_ cxp(-40) 



( 1 + _d$_Y 

\ L T dlna) 



exp(40 Bi3 A' - 40), (9) 



such that the expansion rate is very close to that LCDM at the epoch of BBN where 
the radiation density p r dominates, i.e., no corrections at BBN. Note that TeVeS 
would mimic Dark Matter and Dark Energy if we identify 

P<tF -> PA, (pb + Pr)(r - 1) -> pdm- (10) 
where in the LCDM framework the Hubble expansion is normally modeled with 

3H 2 TT da 

PA + PDM + Pb + Pr = ^-, H=—. (11) 



3, 2008 15:40 WSPC/INSTRUCTION FILE zhaoq2c 



6 HongSheng Zhao 

Using the mirror-imaged a = Lagrangian of the scalar field, we derive the 
following 2nd order ODE for the scalar field <p. 

d 



dt 



/is 



dt 



, o , din a exp(5(/> + d>o)a 3 fdd>\,^. 

= (P6 + „ r )„3, „ = ,„ s E^J (_j (12) 

Note the similarity of this ID equation with the 3D Poisson's equation. 

We can integrate the above equation to solve for (f> as function of In a or the 
physical synchronous time t. We note at current epoch 6 <~ 6o, u s <~ — <~ ^- <~ 

( 2 \-l/2 

(s^r J . where s ~ Vatt^t- 

We then aim to test if the cosmology specified by the above non-fine-tuned 
Lagrangian in TcVeS could matches the behavior of a LCDM universe. We solve 
the equations numerically by iteration of the bare constant G. Assuming a value for 
G, the initial <fi, and are then set by the fact that G e // w G© = 6.67 x 10~ n at 
BBN in order to be consistent with the number of relativistic degrees of freedom at 
temperature of lMeV. The parameters A and K are determined by the boundary 
condition at present day such that we recover the normalisation in the MONDian 
Poisson's equation in galaxies and in solar system. Typically the scalar field tracks 
the matter density, and C and (f> are slow varying functions of redshift. We then 
iterate the parameter G such that the sound horizon angular size at LSS (z = 
1000 — 1100) matches that of LCDM. The parameters typically converge in 20- 
30 iterations. The Hubble constant and cosmic acceleration come out without any 
tuning. 

In Fig. 2a we show a model with a present day matter density Wb = 0.024. 
This is consistent with the baryon density at BBN, hence there is no non-baryonic 
matter in the present model. This model invokes neither cosmological constant nor 
dark matter. The resulting model has Hq = 77km/s/Mpc. The expansion history is 
almost the same as LCDM; slight difference exists in the energy density (hence the 
expansion rate) in the future a > 1. 

To understand whether the above explanation for late acceleration and dark 
matter is unique we have also run models with a more general Lagrangian. 

r< ¥ )=/ ¥ ^ • , » S ^^ T , A = 0. (13) 
Jo 87rG© 1 - as a exp(-0 o ) 

This is so constructed that we recover eq. (1) for any value of a by identifying 
dark matter gravity with the scalar field V<p. This whole sequence of models are 
largely consistent with Dark Matter phenomenology on galaxy scales, with a slight 
preference for a = 1 models in galaxies. Models with non-zero a also have interesting 
effects on the solar system. For example, a model with a = 0.2 would predict (cf. eq. 
1) a constant, non-Keperlian acceleration of ap = a^a^ 1 ~6x 10~ 10 mscc~ 2 in the 
solar system, consistent with the Pioneer Anomaly (although an non-gravitational 
origin is hard to excluded). Such a constant gravity would cause an gravitational 
redshift of 10~ 13 (D/100AU) between the solar system bodies of separation D, 
which could be testable with experiments with accurate clocks in the future (see 
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these proceedings). Calculating the Hubble expansion for models with increasing a, 
we are able to match LCDM in all cases in terms of BBN, LSS, SNela distances, 
and late acceleration. For all these models we have also varied initial conditions 
and found the solutions are stable. The acceleration continues into far future with 
b = acxp((/)) — > est. The amplitude of the modification funtion /i also decreases with 
expansion. Compared to model with a = 0, however, larger a drives up the present- 
day Hubble constant, unless the present day matter density Wb is also increased. 
This is effectively achieved by allowing for relatively massive neutrinos. E.g., for 
a = 0.2 would require the matter density parameter Wb twice the norminal value 

0. 024. implying the need to include massive neutrinos of 0.8 eV. For a = 1 would 
require 2 eV neutrinos as needed for explaining galaxy cluster data (Sanders 2003, 
Pointecouteau & Silk) and the CMB (Skordis ct al. 2006). The latter model is shown 
in Fig.2b. 

5. Conclusion 

In summary, we have focused on one very specific model in the Bckcnstcin theory. We 
have shown that it may be possible to satisfy some of the most stringent cosmological 
observations without the need to introduce/fine-tune dark matter nor dark energy. 
The TeVeS scalar field /z-function (called ^-essence here) can be fixed by galaxy 
rotation curves, and it predicts the right amount of cosmic acceleration, the size 
of the horizon at z = 1000, and the present Hubble constant without fine-tuning. 
The ultimate test of the model should come from simulating the evolution of linear 
perturbations on this background and CMB. By fitting galaxy cluster data and 
the 3rd peak of CMB we could break the degeneracy of models of different a, and 
constrain the neutrino mass. 
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Fig. 1. Left panel shows TeVeS fits to rotation curves of a gas-rich dwarf galaxy NGC1560 and a 
gas-poor larger spiral galaxy NGC4157 (solid curves), adopting ao = 1.2 X 10~ 8 , a = /i-function 
model without neutrinos; the Newtonian rotation curve by baryons for the assumed stellar (M/L)* 
are shown as well (dashed lines). Right panel similar to the left, except for assuming the a = 1 

/^-function and assuming that galaxies are imbeded in a neutrino over-density of 200 X ~ 

2.7 X 10 4 M Q kpc~ 3 or 5000 X ~ 6.7 X 10 5 M Q kpc~ 3 (the two values brackets the typical gas 
density of x-ray clusters on average and in the centres). The Newtonian rotation curves of the 
constant neutrino cores are also shown (dotted lines). 
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Fig. 2. Shows the values for the TeVeS ao parameter derived to fit individual strong lensing 
Einstein radius of 50 CASTELS multi-imaged systems, assuming M*/L* = 4 (circles); also shown 
are the effects of raising/lowering M*/L* by a factor of 2 (solid vertical lines) or a factor of 4 
(dotted vertical lines). A few outliers are labeled. The right panel shows the Newtonian acceleration 
GM„/R 2 E vs. the critical gravity (related to the critical surface density (c 2 /4nGDi)(D a / Di a ) in 
GR) is the mininal local gravity for a lens to form Einstein rings. The dashed line is a prediction 
for point lenses in TeVeS a = model. 
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TeVeS (solid) mimics LCDM (dashed) TeVeS (solid) mimics LCDM (dashed) 

i 1 t — ' ' ' w i — — — — a in' t — ' ' \n — — — — a 




Scale factor a Scale factor a 



Fig. 3. compares ACDM (dashed) with Zero-A-TcVcS flat cosmologies (solid) (left panel assumes 
zero mass for neutrinos and a ^t-essence with a = 0; right panel assumes 2eV neutrinos and a = 1 
model). Shown are the co-moving distance D com vs. the physical scale factor a in log-log diagram 
ovcrplotted with SNIa data (small symbols) up to redshift 2. Likewise shows the horizon, the 
Hubble parameter H in units of (Mpc _1 c) in two theories. The evolution of the scalar field <f> 
and /i can be inferred from (thin solid lines) aexp(</>) and with the cutoff of = 0.005 be 
adopted for numerical reasons. 



